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ABSTRACT 

We develop the theory of jitter radiation from GRB shocks containing small-scale magnetic fields 
and propagating at an angle with respect to the line of sight. We demonstrate that the spectra vary 
considerably: the low-energy photon index, a, ranges from to —1 as the apparent viewing angle goes 
from to 7r/2. Thus, we interpret the hard-to-soft evolution and the correlation of a with the photon 
flux observed in GRBs as a combined effect of temporal variation of the viewing angle and relativistic 
aberration of an individual thin, instantaneously illuminated shell. The model predicts that about a 
quarter of time-resolved spectra should have hard spectra, violating the synchrotron a = —2/3 line of 
death. The model also naturally explains why the peak of the distribution of a is at a « —1. The 
presence of a low-energy break in the jitter spectrum at oblique angles also explains the appearance of 
a soft X-ray component in some GRBs and a relatively small number of them. We emphasize that our 
theory is based solely on the first principles and contains no ad hoc (phenomenological) assumptions. 

Subject headings: gamma rays: bursts — radiation processes — shock waves — magnetic fields 



1. INTRODUCTION 

Rapid spectral variability is a remarkable, yet unex- 
plained feature of the prompt GRB emission. The vari- 
ation of the hardness of the spectrum and the hard- 
to-soft evolution are the most acknowledged features 
fBhat, ct al.' "1994"; 'Crider. et all Il997t iFrontera. erall 
.2000: Rvde fc Pc trosian 2003i). ^ quite remarkable "track- 
ing" behavior, when the low-energy spectra l index a fol- 
lows (or correlates with) the photon flux ijCrider. et alJ 
Il997|) is particularly intriguing. An e xample of such tren d 
is shown in Fig. ^(we used data from IPreece. et al]l2000j) . 

In this paper, we demonstrate that the spectral index- 
flux correlation (or the hardness-intensity correlation) is 
a natural and inevitable prediction of the coUisionless 
relativistic shock model and the jitter radiation mech- 
anism from small-scale s hock-ge nerated magnetic fields 
ijMedvedev fc Loeblll999t lMedvedev.2000 ). Quite possi- 
bly, th e appearance of a soft X-ray component in some 
GRBs dp reece. can also be interpreted within 

this model. The theory yields a right number of the 
synchrotron- violating GR Bs, i.e., with a > —2/3, ijKatd 
11994 IPreece. et al.lll998l) and natura lly explains why th e 
majority of GRBs have a around —1 l|Preece. et aD|2Q03). 
We emphasize that our theory is based solely on the first 
principles and contains no ad hoc assumptions. 
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where we neglect plasma dispersion, which is vanishing 
for the GRB photons. Here 7 is the Lorentz factor of 
a radiating particle and w^j' = J we'" * dt is the Fourier 
component of the transverse particle's acceleration due to 
the Lorentz force, F^. This temporal Fourier transform is 
taken along the particle trajectory, r = Tq -I- vi. It should 
appropriately be expressed via the statistical properties of 
the magnetic field. To make this sect ion self-cont ained, we 
will, in part, follow the derivation of iFleishmanI |2005j). 

We need to express the temporal Fourier component of 
the acceleration w = FlI^vti taken along the particle tra- 
jectory in terms of the Fourier component of the field in the 
spatial and temporal domains. Taking the Fourier trans- 
form of w(ro -l- vi, f), we have 



= (2^)- 
= (2^)- 



e*'^'* dt 



wo,k<5(w' -n + k-^^) e'^-"" dndk, (2) 



2. THEORY OF JITTER RADIATION IN 3D 

The angle-averaged spectral power emitted by a rela- 
tivistic particle moving through small-scale random mag- 
netic fields, under the assumption that the defiection angle 
is negligible and the particle trajectory i s a straight line, 
has been derived elsewhere ( Rybicki & LightmanI 119791: 
ILandau fc Lifshitdll97ll: lMedvedev..2OO0^) and it reads 

'^also at the Institute for Nuclear Fusion, RRC "Kurchatov Institute" 
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where we used that / e*"* = 2Tr5{uj). In a statisti- 
cally homogeneous turbulence, |w^/p should not depend 
on the initial point, Tq, of the particle trajectory. There- 
fore, we average it over rg as (|w^'p) = J |w^'p(iro, 
where V is the volume of the spatial domain. Using that 
/ e*(k-ki) ro ^j.^ ^ (27r)3,5(k - ki), we finally have: 



|wn,kp'5(c^' -n + k-v) dndk. 



(3) 

The Lorentz acceleration, w = {e/jmc)^ x B, can be 
written as Wa = {e/jmc)^ea0j{vpBj — v^Bfi). Using the 
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identity, eaf3^ea\^, = SpxS^t^ - Spf^S^x, we obtain 

|wa,kP = {ev/jmc)^{So.p - v~^Vo.v„) i?8,kSo^- (4) 

In a statistically homogeneous random magnetic field, the 
tensor k^ok '^^^ expressed via the Fourier trans- 
form of the field correlation tensor 

Bn^Bn^ = TVK^p{n,k) = Tvje'^'''-'^-^^ Ko.i3{r,t) drdt, 

where T is the size of the temporal domain and Kapir, t) — 
T-iy-i/ Ba{v' ,t')Bp{Y' + Y,t' + t)dY'dt is t he second- 
order correlation tensor of the magnetic field l)Fleishmanl 
120051) . 

In the static case, i.e., when the magnetic field is inde- 
pendent of time, Eqs. l|2Jl, Q read as 



(|w^f ) = (2^F)-i j + k • v) dk. 



(6) 



|wkP = {ev/^mcY{6ap - V^'^VaVp) VKap{\<i). (7) 
3. THE MAGNETIC FIELD SPECTRUM 

We adopt the following geometry: a shock is lo- 
cated in the cc-y-plane and is propagating along z- 
direction. As it has in itially been demonstrated by 
iMedvedev &"Loe 9 (TOQ?) and later confirmed via 3D 
PIC simulations (Silva, et al. 2003; Nishikawa, ct al. 2003; 
iFrederiksen. et al1l2004() . the magnetic field at relativis- 
tic shocks is described by a random vector field in the 
shock plane, i.e., the x-y-plane. As the shock is prop- 
agating through a medium, the produced field is trans- 
ported downstream (in the shock frame) whereas new field 
is continuously generated at the shock front. Thus, the 
field is also random in the parallel direction, i.e., the z- 
direction. Thus, Weibel turbulence at the shocks is highly 
anisotropic. Both the theoretical considerations and re- 
alistic 3D simulations of relativistic shocks indicate that 
the dynamics of of the Weibel magnetic fields in the shock 
plane and along the normal to it are decoupled. Hence, 
the Fourier spectra of the field in the x — y plane and in z 
direction are independent. Thus, for the Weibel fields at 
shocks, the correlation tensor has the form 



(8) 



where n is the unit vector normal to the shock front, C is 
the normalization constant proportional to (i?^), fz and 
fxy are the magnetic field spectra along n and in the shock 
plane, respectively, kj_ = (fc^ -I- A:^)^/^ and /c|| = k^, and fi- 
nally, the tensor {Sap — naTip) is symmetric and its product 
with n is zero, implying ort hogonality of n and B. 

Numerical simulatrions ijFrederiksen. et all l2004j) also 
indicate that the field transverse spectrum, fxy^ is well de- 
scribed by a broken power-law with the break scale compa- 
rable to the skin depth, c/wp, where ujp — {ine^n/Tmy^^ 
is the relativistic plasma frequency and F is the shock 
Lorentz factor. We expect that the spectrum f^, has sim- 
ilar properties. Therefore, we use the following models: 



fxy{k±) 



k 



(9) 



where K|| and are parameters (bei ng, in general, a 
funct ion of the distance from the front, IMedvedev. et alJ 
120051) determining the location of the peaks in the spec- 
tra, tti, 02, /32 are power-law exponents below and 
above a spectral peak (/3i > a2 + 1/2 and /32 > ^2 + 1, for 
convergence at high-fc). Note that /3i^2 — > oo corresponds 
to spectra with a sharp cut-off. The asymptotes of these 
functions are 



/(fc)a 



if A; <i; K, 
2/3, if fc > K. 



(10) 



4. RADIATION SPECTRA FROM A SHOCK VIEWED AT 
DIFFERENT ANGLES 

We now evaluate Eqs. ©,10 ■ The scalar product of 
the two tensors is 

{^aS3-VaVl3lv^){6aS3-nanii) = X^inaV^Y jv^ = l+COS^ 6, 

(11) 

where we used that 6aa = 3. Here 9 is the angle be- 
tween the normal to the shock and the particle velocity (in 
an observer's frame), which is approximately the direction 
toward an observer, that is v[|k for an ultra-relativistic 
particle (because of relativistic beaming, the emitted ra- 
diation is localized within a narrow cone of angle ~ 1/7)- 
Eq. ^ becomes 

(|w^f ) = ^(l+COS^e) jUk^^)Uy{k^)S{uj' + -k-v) dk\\d^ki 

(12) 

Equations lfT|l. (|12() fully determine the spectrum of jitter 
radiation from a GRB shock. We now consider special 
cases. 

4.1. A shock viewed head-on, = 

For a shock moving towards an observer, n|jk, hence 
8 = (because n||v and k||v for 7^1) and k • v = kzV 
(of course, u « c). Equation (|12|l then reads 



{\^u.'?)^^-^C fxy / fz{kz)5{Lu' + kzV)dk 



iCMv\)fxyM~L0'/v), 



(13) 



where fxy = / fxy{k±)d^k±. A pparently, t his ca se is anal- 
ogous to the ID model used bv IMedvedev! 1)20001) . We now 
determine low-cj and high-w asymptotics for the spectrum 
given by Eq. @. In order to simplify the analysis, we 
neglect the small second and third terms in the brackets 
in Eq.JQ and assume that ai > 1/2, /3i > 1. This slightly 
changes the shape of the radiation spectrum near a peak, 
but does not affect the asymptotic behavior. We have 



dW 

duj 



OC LjJ 



2ai-2/3i 



y 



2ai 



dy 



(14) 



where y — uj' /uj, ??|| = k^v/lo, and ya = 1/27^. 

At low frequencies, 77|| ^ yo, (that is, uj ^ k^^vj'^), 
the right hand side of Eq. (|14|) can approximately 
be evaluated as RHS - w2ai-2ft(jj;ii ^201-2^-2/31^^ _|_ 

OC At 



UJ 

/^^2/2ai-2-2/3i^^) OC ^2ai-2/3i„2ai-l-2/3i _ , ,i 



'»7|| 



high frequencies 7711 <^ yo, (that is, w ^ K||-y7 ), the 
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right hand side becomes (note, yo — const.) RHS ^ 
^2ai-2/3i j°Oy2ai-2~2/3i^j^ ^ ^2ai-2/3i^ Combining the 



'ya 

results, we conclude that 



dW 
duj 



OJ 



,2ai-2/3i 



if < K||U7^, 

if ^ K||U7^. 



(15) 



4.2. ^ shock viewed edge-on, Q =^ Tr/2 

An ultra-relativistic shock moving at an angle ~ l/F 
with respect to the line of sight is seen nearly edge-on be- 
cause of relativistic aberration. In this case, the shock is 
seen as if n ± k. Therefore Q = tt/2 and k ■ v = kxV, 
where we assumed that an observer is located on the x- 
axis. Equation H12|) then becomes 

(|w^,|'> = {2Tr)-'CTz f Uyik±)diu' + k^v) dk^dky 



(C/27r|«|)/, / f,y{^iLo'/vr + kl)dky 
(l+y2)a. 



oc (w') 



/\2a2-2/32 + l 



ivi + '^ + y 



dy, (16) 



where fz = ! fz{k\\)dk\\. In the last line, we introduced 
y = kyv/uj' and ry* = k±_v/uj' . The integral in Eq. (|16|l 
is independent of w' when ry* -C 1, i.e., at high frequen- 
cies, u' » K±v. At low frequencies, ?/* ^ 1, the inte- 
gral is dominated by y ~ 77* (as in Eq. d]), hence it is 



oc r;^"=^-2/32+i ^ ^/^(2a2-2/32+i), Thus, wc havc 



(^') 
(c') 



/\2a2-2/32 + l 



if <C K±f , 

if 3> Kj^f . 



(17) 



Comparing Eqs. ((T7|l and ((Tn|l . we approximate (|wtj/p) 
by a function as in Eq. © with a = and —2/3 = 
2a2 — 2/32 + 1- We can now find the asymptotes of the 
spectra from Eq. (Q) . The analysis is analogous to that of 
Eq. (dl) and it yields 



dW 
duj 



e=7r/2 



,202-2/32 + 1^ ifcc;>K^W72. ^^^> 



This result is analogous to that o f iFleishmanI HHoH). 

4.3. A shock viewed at oblique angles, < < 7r/2 

For a shock viewed at an oblique angle, k = x/c sin Q 
zk cos Q, we have k • v = k^v sin Q + k^v cos Q. Hence 



(|w.f) 



2^ 



fz{kz)fxyikx)Siuj' + k • v) dkzdk^ 
/z 1 "~7T + tane^ fxy{kx) dk^, 



cx (w') 



cos 

(l + 2/)2"i 



N? + (i + y)Y^ {ril + y^Y^ 



(19) 



where we used that f^^—kz) = fz{kz) and we defined C* = 
C/(2^|wcose|)and/,y(fc,) = //,j,((fc2+fc2)i/2)^fc^^ One 

can show that this function is very similar to fxy{k\_) with 
fc_L replaced with k^^ so we use it in our analysis. In the 
last hne we introduced C, = 2ai — 2/3i -I- 2a2 — 2/32 + 1, 



y = kxVsmQ/uj' , rji — k\\v cos / lu' , 772 = K±vsmQ/uj'. 
If K±_ ^ K|| and 0^1, then 771 ^ 772 and we can con- 
clude that (|wtj/p) has two breaks: uj[ = KyticosG and 
^'2 = K±vsiTiQ. A more detailed numerical analysis and 
the resulting radiation spectrum are discussed below. 

5. INTERPRETATION OF PROMPT GRB SPECTRA 

In the standard internal shock model, each emission 
episode is associated with illumination of a thin shell, 
— an internal shock and the hot and magnetized post- 
shock material. We assume that the shell is spherical (at 
least within a cone of opening angle of ^ l/F around 
the line of sight) and this shell is simultaneously illu- 
minated for a short period of time. The observed pho- 
ton pulse is broadened because the photons emitted from 
the patches of the shell located at larger angles, 1} . from 
the li ne of sight arrive at progressively later times l)Piranl 
Il999j) . This effect naturally explains the fast-rise-slow- 
decay lightcurves of individ ual pulses l|Rvde fc PetrosianI 
l2002t iKocevski. et al.ll2003l) . Because of relativistic aber- 
ration, the apparent viewing angle, Q, is greater than 
and approaches 8 ~ 7r/2 (the shell is seen edge-on) when 
d ~ l/F. Thus, there must be a tight correlation between 
the observed spectrum and the observed photon flux, be- 
cause they are, in essence, different manifestations of the 
same relativistic kinematics effect. 

Let us now discuss specific properties of the predicted 
spectra. Fig. El represents full numerical solutions of Eqs. 
Ql, ®i l|12|l for three different viewing angles. In cal- 
culation of dW/dio, the emitting electrons were assumed 
monoenergetic, for simplicity. The extension to a stan- 
dard power-law with a sharp low-energy cutoff, N ^ 7"^ 
for 7 > 7niin is straightforward: the low-o; spectral slope 
remains unchanged, and the high-o; slope is equal to ^ 
or — (p — l)/2, whichever is greater (neglecting e~ cool- 
ing). An important fact to note is that the jitter radiation 
spectrum varies with the viewing angle. When a shock ve- 
locity is along the line of sight, the low-energy spectrum is 
hard Fi, (x i^^, harder than the "synchrotron line of death" 
(Fjy oc 7/^/^). As the viewing angle increases, the spectrum 
softens, and when the shock velocity is orthogonal to the 
line of sight, it becomes oc 7/°. Another interesting fea- 
ture is that at oblique angles, the spectrum does not soften 
simultaneously at all frequencies. Instead, there appears a 
smooth spectral break, which position depends on O. The 
spectrum approaches ~ below the break and is harder 
above it. This softening of the spectrum at low i^'s could 
be interpreted as the appearance of an additional soft X- 
ray component, sim ilar to that found in some of GRBs 
ijPreece. et al.lll995|) . 

Fig. 131 represents the spectral slope evaluated at fre- 
quencies about 10 and 30 times below the spectral peak. 
These frequencies correspond to the edge of the BATSE 
window for bursts with the peak energy of about 200 keV 
and 600 keV, respectively. Hence, the spectral slope, 
acRB, will be close to those obtained from the data fits. 
Since Q{t) increases with time during an individual emis- 
sion episode, the curves roughly represent the temporal 
evolution of acRB- Assuming that time-resolved spec- 
tra are homogeneously distributed over 6, one can es- 
timate the relative fraction of the synchrotron-violating 
GRBs (i.e., those with acRB + 1 > 1/3) as about 25%, 
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which is very close to the 30% obtained from the data 
(jPreece. et al.ll2000l) . Most of the GRBs, -75%, should, 
by the same token, be distributed around acRB ~ — 1- 
Note also that time-integrated GRB spectra should have 
Q^GRB around minus one, as well. This explains why 
a relatively large sample of synchrotron-violating bursts 
is present in the tinie-resolved BSAX and BATSE data 
l|Preece. et al.ll2000t lEYontera. et al]l2000tl . Since time- 
integrated data are dominated by a ~ —1 spectra, the 
synchrotron-violating GRBs should be practically absent 
from the time-intergtated i3^ Tgi? and HETE-II spectral 
data ijBarraud. et al.ll2003(l . We stress that the question 
of why the peak of the acRB-distribution is at q;grb — —1 
and not at some other "physically motivated" value of 
—3/2 or —2/3, has had no satisfactory explanation un- 
til now. Finally, the spectral softening which looks like 



an additional soft X-ray component should appear within 
the detector spectral window when 9 20°± few degrees. 
Thus, we estimate that this X-ray excess can be detectable 
in about 10% o f GRBs, which is ag ain quite close to the 
observed 15% (jPreece. et al.lll995j) . Of course, a care- 
ful statistical analysis, which takes into account uneven 
sampling (more time-resolved spectra for brighter parts of 
the bursts), statistical and systematic errors, biases intro- 
duced by fits to particular spectral models, etc., is very de- 
sirable. However, the very fact that relative sizes of GRB 
populations fall in the right bulk part is very encouraging. 

This work has been supported by NASA grant NNG- 
04GM41G, DoE grant DE-FG02-04ER54790, and the KU 
GRF fund. 
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Fig. 1. — The normalized photon flux (diamonds) and the low-energy power-law index a (squares) vs time for BATSE trigger #2953. The 
data are from the time-resolved spectral fits bv iPreccc. et al] J2000i) . 
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Fig. 2. — The log — log plots of |wkP vs k {thin lines) and of dW/doj vs ui (thick lines), for three viewing angles © = 0, tt/IO, it/2. The axes 
units are arbitrary. In this calculation we used fz = f^y with a = 2, 13 = 20, k = 1, v = 1. The exponent = f{o,/3) is model dependent 
[c.f., Eqs. I15i . 1181 1. We also chose 7 = 1 in order to align the peaks of |wkP and dW/duj. Note that the actual peaks are at values k, u) 
lower than 10 by a factor two or three. Note also that the spectrum dW/du) levels off at oblique angles at frequencies much smaller than 
Ki;7^ sinS, whereas |wk|^ indeed starts to flatten at A; ~ /ctisinS. 
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Fig. 3. — The low-energy spectral slope evaluated at oj ^ 0.03 [solid curve) and oj 0.01 [dashed curve) in units of Fig. |3 i.e., at oj's 
about 10 and 30 times below the spectral peak, corresponding to the \aw-u) edge of the BATSE window for GRBs with Bp of about 200 and 
600 keV, respectively. Note that abour 25% of spectra violate the synchrotron line of death (solid horizontal line). 



